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Abstract 

We investigate a functional limit theorem (homogenization) for Reflected Stochastic Dif- 
ferential Equations on a half-plane with stationary coefficients when it is necessary to analyze 
both the effective Brownian motion and the effective local time. We prove that the limiting 
process is a reflected non-standard Brownian motion. Beyond the result, this problem is 
known as a prototype of non-translation invariant problem making the usual method of the 
"environment as seen from the particle" inefficient. 
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1 Introduction 

Statement of the problem 

This paper is concerned with homogenization of Reflected Stochastic Differential Equations 
(RSDE for short) evolving in a random medium, that is (see e.g. [T2] ) 

Definition 1.1. (Random medium). Let (Q,Q,(i) be a probability space and {r x ;x £ R d } be a 
group of measure preserving transformations acting ergodically on 57, that is: 

1) VA G G,Vx £ R d , h{t x A) = n{A), 

2) If for any x £ R d , t x A = A then fj,(A) = or 1, 

3) For any measurable function g on (Q, Q, (J,), the function {x, ui) <— > g(r x co) is measurable on 

(w l x n,B(R d )®g). 

The expectation with respect to the random medium is denoted by M. In what follows we 
shall use the bold type to denote a random function g from fl x W into M. n (n > 1 and p > 0). 

A random medium is a mathematical tool to define stationary random functions. Indeed, 
given a function / : il — > K, we can consider for each fixed u> the function x G K d i— > f(r x uj). This 
is a random function (the parameter uj stands for the randomness) and because of 1) of Definition 
11.11 the law of that function is invariant under K d -translations, that is both functions f(r.ui) and 
f(T y +.u>) have the same law for any y G M. d . For that reason, the random function is said to be 
stationary. 

We suppose that we are given a random d x <i-matrix valued function er : £1 — > R dxd , two 
random vector valued functions b, 7 : Q —> and a d-dimensional Brownian motion B defined 
on a complete probability space (Q!,!F, P) (the Brownian motion and the random medium are 



independant). We shall describe the limit in law, as e goes to 0, of the following RSDE with 
stationary random coefficients 

(1) dX £ t = e" 1 b(r XtVe w) dt + (t(t x , /e uj) dB t + 70xf/ e w) dK e tl 

where X e ,K £ are (J-'t) i-adapted processes (JFt is the cr-field generated by B up to time t) with 
constraint Xf £ D, where D C M. d is the half-plane {(xi, . . . , Xd) £ K d ; x\ > 0}, K e is the so-called 
local time of the process X e , namely a continuous nondecreasing process, which only increases on 
the set {t; Xf £ dD}. The reader is referred to [14] for strong existence and uniqueness results to 
|T]) (see e.g [23] for the weak existence), in particular under the assumptions on the coefficients <t, b 
and 7 listed below. Those stochastic processes are involved in the probabilistic representation of 
second order partial differential equations in half-space with Neumann boundary conditions (see 
|18] for an insight of the topic). In particular, we are interested in homogenization problems for 
which it is necessary to identify both the homogenized equation and the homogenized boundary 
conditions. 

Without the reflection term j(X^/e) dK%, the issue of determining the limit in (fT]) is the 
subject of an extensive literature in the case when the coefficients b, a are periodic, quasi-periodic 
and, more recently, evolving in a stationary ergodic random medium. Quoting all references is 
beyond the scope of this paper. Concerning homogenization of RSDEs, there are only a few works 
dealing with periodic coefficients (see [TJ [2j [3] |22] ) . As pointed out in [2], homogenizing (fT]) in a 
random medium is a well-known problem that remains unsolved yet. There are several difficulties 
in this framework that make the classical machinery of diffusions in random media (i.e. without 
reflection) fall short of determining the limit in (fT]). In particular, the reflection term breaks the 
stationarity properties of the process X e so that the method of the environment as seen from 
the particle (see [TB] for an insight of the topic) is inefficient. Moreover, the lack of compactness 
of a random medium prevents from using compactness methods. The main resulting difficulties 
are the lack of invariant probability measure (IPM for short) associated to the process X s and 
the study of the boundary ergodic problems. The aim of this paper is precisely to investigate the 
random case and prove the convergence of the process X 6 towards a reflected Brownian motion. 
The convergence is established in probability with respect to the random medium and the starting 
point x. 

We should also point out that the problem of determining the limit in (fT]) could be expressed 
in terms of reflected random walks in random environment, and remains quite open as well. In 
that case, the problem could be stated as follows: suppose we are given, for each z £7L d satisfying 
\z\ = 1, a random variable c(-, z) : f2 — >]0; +oo[. Define the continuous time process X with values 
in the half-lattice L = N x Z d_1 as the random walk that, when arriving at a site x £ L, waits 
a random exponential time of parameter 1 and then performs a jump to the neighboring sites 
y £ L with jump rate c(T x u>,y — x). Does the rescaled random walk eX t j e i converge in law 
towards a reflected Brownian motion? Though we don't treat explicitly that case, our proofs can 
be adapted to that framework. 

Structure of the coefficients 

Notations: Throughout the paper, we use the convention of summation over repeated indices 
Si=i °i-di = c-idi and we use the superscript * to denote the transpose A* of some given matrix 
A. If a random function (p : — ► K possesses smooth trajectories, i.e. for any uj £ ft the mapping 
x £ M. d i— > (f(T x uj) is smooth with bounded derivatives, we can consider its partial derivatives at 
denoted by Di(p, that is Di<p(uj) — d Xi (x t— > (p(T x u))i x=Q . □ 

We define a = crcr* . For the sake of simplicity, we assume that £ £1 the mapping x £ 
M. d \— > <j(t x lq) is bounded and smooth with bounded derivatives of all orders. We further impose 
these bounds do not depend on lu. 

Now we motivate the structure we impose on the coefficients b and ~y. A specific point in 
the literature of diffusions in random media is that the lack of compactness of a random medium 
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makes it impossible to find an IPM for the involved diffusion process. There is a simple argument 
to understand why: since the coefficients of the SDE driving the Revalued diffusion process are 
stationary, any R d -supported invariant measure must be stationary. So, unless it is trivial, it 
cannot have finite mass. That difficulty has been overcome by introducing the "environment 
as seen from the particle" (ESFP for short). It is a il- valued Markov process describing the 
configurations of the environment visited by the diffusion process: briefly, if you denote by X the 
diffusion process then the ESFP should match tx^j- There is a well known formal ansatz that 
says: if we can find a bounded function / : Q — > [0, +00 [ such that, for each u) € fi, the measure 
f{T x uj)dx is invariant for the diffusion process, then the probability measure f(u>)dfx (up to a 
renormalization constant) is invariant for the ESFP. So we can switch an invariant measure with 
infinite mass asociated to the diffusion process for an IPM associated to the ESFP. 

The remaining problem is to find an invariant measure (of the type f(r x ui)dx) for the diffusion 
process. Generally speaking, there is no way to find it excepted when it is explicitly known. In 
the stationary case (without reflection), the most general situation when it is explicitly known is 
when the generator of the rescaled diffusion process can be rewritten in divergence form as 

(2) £ s f = \e 2V ^'^d Xi {e- 2V ^/^\ aij +H ij ){T x/e u>)d x j), 

where V : £1 — * R is a bounded scalar function and H : £1 — * W lxd is a function taking values in 
the set of antisymmetric matrices. The invariant measure is then given by e^^/^'dx and the 
IPM for the ESFP matches e 2V ^dfi. However, it is common to assume V = H = to simplify 
the problem since the general case is in essence very close to that situation. Why is the existence 
of an IPM so important? Because it entirely determines the asymptotic behaviour of the diffusion 
process via ergodic theorems. The ESFP is therefore a central point in the literature of diffusions 
in random media. 

The case of RSDE in random media does not derogate this rule and we are bound to find 
a framework where the invariant measure is (at least formally) explicitly known. So we assume 
that the entries of the coefficients b and 7, defined on CI, are given by 

(3) Vj = 1, . . . , d, b 3 = -Didij, = dji. 

With this definition, the generator of the Markov process X e can be rewritten in divergence form 
as (for a sufficiently smooth function / on D) 

(4) C £ f = -d Xi (a t] {T x/e uj)d x J) 

with boundary condition ~fi(T x / E ui)d Xi f — on dD. If the environment u> is fixed, it is a simple 
exercise to check that the Lebesgue measure is formally invariant for the process X e . If the ESFP 
exists, the aforementioned ansatz tells us that /i should be an IPM for the ESFP. Unfortunately, 
we shall see that there is no way of defining properly the ESFP. The previous formal discussion 
is however helpful to provide a good intuition of the situation and to figure out what the correct 
framework must be. Furthermore the framework (|3|) also comes from physical motivations. As 
defined above, the reflection term 7 coincides with the so-called conormal field and the associated 
PDE problem is said to be of Neumann type. From the physical point of view, the conormal 
field is the "canonical" assumption that makes valid the mass conservation law since the rela- 
tion aji(r x f e uj)d Xj f = on dD means that the flux through the boundary must vanish. Our 
framework for RSDE is therefore to be seen as a natural generalization of the classical stationary 
framework. 

Remark. It is straightforward to adapt our proofs to treat the more general situation when the 
generator of the RSDE inside D coincides with ^ . In that case, the reflection term is given by 

I3 = "> 11 
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Without loss of generality, we assume that an = 1. We further assume that a is uniformly 
elliptic, i.e. there exists a constant A > such that 

(5) efl, AI < a(u) < A- 1 !. 

That assumption means that the process X e diffuses enough, at each point of D, in all directions. 
It is thus is a convenient assumption to ensure the ergodic properties of the model. The reader 
is referred, for instance, to [5J [2U1 HI] f° r various situations going beyond that assumption. We 
also point out that, in the context of RSDE, the problem of homogenizing (TTJ) without assuming 
([5j becomes quite challenging, especially when dealing with the boundary phenomena. 

Main Result 

In what follows, we indicate by PJ the law of the process X £ starting from x E D (keep in mind 
that this probability measure also depends on to though it does not appear through the notations). 
Let us consider a nonnegative function \ : D — > R_|_ such that Jq xi x ) dx = 1. Such a function 
defines a probability measure on D denoted by x(dx) = x(x)dx. We fix T > 0. Let C denote the 
space of continuous D x R + -valued functions on [0, T] equipped with the sup-norm topology. We 
are now in position to state the main result of the paper: 

Theorem 1.2. The C-valued process (X e ,K e ) £ converges weakly, in H®X probability, towards 
the solution (X, K) of the RSDE 

(6) X t = x + A^ 2 B t + TK t , 

with constraints X t € D and K is the local time associated to X . In other words, for each bounded 
continuous function F on C and 5 > 0, we have 

lim (i® x{(u,x) e Q. x D; \E S JF(X £ , K e )) - E X (F(A > , K))\ > 5} = 0. 

The so-called homogenized (or effective) coefficients A and T are constant. Moreover A is invert- 
ible, obeys a variational formula (see subsection \2.5\ for the meaning of the various terms) 

A= inf M[(I + Dip)*a(l + Dip)] , 

and r is the conormal field associated to A, that is Ti = An for i = 1, . . . , d. 

Remark and open problem. The reader may wonder whether it may be simpler to consider the 
case 7,- = Su where 6 stands for the Kroenecker symbol. In that case, 7 coincides with the normal 
to 3D. Actually, this situation is much more complicated since one can easily be convinced that 
there is no obvious invariant measure associated to X e . 

On the other side, one may wonder if, given the form of the generator inside D, one 
can find a larger class of reflection coefficients 7 for which the homogenization procedure can be 
carried through. Actually, a computation based on the Green formula shows that it is possible 
to consider a bounded antisymmetric matrix valued function A : — > M. dxd such that Ay = 
whenever i = 1 or j = 1, and to set ~fj = aj± + DiAji. In that case, the Lebesgue measure 
is invariant for X e . Furthermore, the associated Dirichlet form (see subsection \2.3\) satisfies a 
strong sector condition in such a way that the construction of the correctors is possible. However, 
it is not clear whether the localization technique based on the Girsanov transform (see Section \2.1\ 
below) works. So we leave that situation as an open problem. 

The non-stationarity of the problem makes the proofs technical. So we have divided the 
remaining part of the paper into two parts. In order to have a global understanding of the proof 
of Theorem 11.21 we set the main steps out in Section [2] and gather most of the technical proofs 
in the Appendix. 
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2 Guideline of the proof 



As explained in introduction, what makes the problem of homogenizing RSDE in random medium 
known as a difficult problem is the lack of stationarity of the model. The first resulting difficulty 
is that you cannot define properly the ESFP (or a somewhat similar process) because you cannot 
prove that it is a Markov process. Such a process is important since its IPM encodes what the 
asymptotic behaviour of the process should be. The reason why the ESFP is not a Markov process 
is the following. Roughly speaking, it stands for an observer sitting on a particle X£ and looking 
at the environment tx'OJ around the particle. For this process to be Markovian, the observer 
must be able to determine, at a given time t, the future evolution of the particle with the only 
knowledge of the environment tx^ui. In the case of RSDE, whenever the observer sitting on the 
particle wants to determine the future evolution of the particle, the knowledge of the environment 
Tjsqw is not sufficient. He also needs to know whether the particle is located on the boundary dD 
to determine if the pushing of the local time K\ will affect the trajectory of the particle. So we 
are left with the problem of dealing with a process X e possessing no IPM. 

2.1 Localization 

To overcome the above difficulty, we shall use a localization technique. Since the process X e is 
not convenient to work with, the main idea is to compare X e with a better process that possesses, 
at least locally, a similar asymptotic behaviour. To be better, it must have an explicitly known 
IPM. There is a simple way to find such a process: we plug a smooth and deterministic potential 
V : D — > M into ([4| and define a new operator acting on C 2 (D) 

p 2V(x) d 

(7) C v = J2 d Xz {e' 2V ^ ai] {T x/e uj)d X] ) = - d Xl V(x)a tJ (T x/£ cj)d Xj , 

with the same boundary condition ■y i (r x / E uj)d Xi — on dD. If we impose the condition 

(8) / e 2V ^dx = l 

J D 

and fix the environment to, we shall prove that the RSDE with generator C v inside D and 
boundary condition ~fi{T x / e to)d Xi = on dD admits e 2V ^dx as IPM. 

Then we want to find a connection between the process X e and the Markov process with 
generator C v inside D and boundary condition ■~f i {T x / e uj)d Xi — on dD. To that purpose, we 
use the Girsanov transform. More precisely, we fix T > and impose 

(9) V is smooth and d x V is bounded. 

Then we define the following probability measure on the filtered space (CI'; T, (^ r t)o<t<T) 
d¥%* = exp ( - J d Xi V{X e r )* i:j {Txs r / e u) dBJ r-\J Q ^V{Xf.)a ij {T X . /e u))d 3lj V{X%) drj dF%. 

Under ¥%* , the process £? t * = B t + f Q a-(T X ^/ £ uj)d x V(X^) dr (0 < t < T) is a Brownian motion 
and the process X e solves the RSDE 

(10) dXl = E- 1 6(T JCf/e w) dt - a{T XUe uj)d x V(XI) dt + v(T Xt/e u) dB* t + 7 ( Tjf . /e w) dKf 

starting from X$ — x, where K £ is the local time of X £ . It is straightforward to check that, 
if B* is a Brownian motion, the generator associated to the above RSDE coincides with ([7]) 
for sufficiently smooth functions. To sum up, with the help of the Girsanov transform, we can 
compare the law of the process X e with that of the RSDE (fT0|) associated to C v . 



5 



We shall see that most of the necessary estimates to homogenize the process X e are valid under 
PJ*. We want to make sure that they remain valid under P|. To that purpose, the probability 
measure P| must be dominated by P|* uniformly with respect to e. From ([9]), it is readily seen 

that C = sup e>0 (E e x *[(-§kr) 2 ]) 1/2 < +00 (C only depends on T,\a\oo and sup^, \d x V\). Then 
the Cauchy-Schwarz inequality yields 

(11) Ve > 0, VA .^-measurable subset , f%{A) < C(W%*{A)) 1/2 . 

In conclusion, we summarize our strategy: first we shall prove that the process X e possesses 
an IPM under the modified law P e *, then we establish under P e * all the necessary estimates to 
homogenize X e , and finally we shall deduce that the estimates remain valid under P e thanks to 
(fTTj) . Once that is done, we shall be in position to homogenize ((T|). 

To fix the ideas and to see that the class of functions V satisfying ([8|) © is not empty, we 
can choose V to be equal to 

(12) V(xi, . . . , x d ) = Axi + A(l + x\ + ■ ■ ■ + x 2 d ) 1/2 + c, 

for some renormalization constant c such that fg e~ 2V ( X ^ dx = 1 and some positive constant A. 

Notations for measures. In what follows, P e (resp. P e *) stands for the averaged (or annealed) 
probability measure M J 5 ¥%{-)e- 2V ^ dx (resp. M f D W*{-)e~ 2V ^ dx), and W (resp. W*) for 
the corresponding expectation. 

¥* D and Pg D respectively denote the probability measure e~ 2VtyX ^ dx ® dfi on D x VI and the finite 
measure e~ 2V ^ dx ® d/j, on dD x Q. M* D and M.g D stand for the respective expectations. 

2.2 Invariant probability measure 

As explained above, the main advantage of considering the process X e under the modified law 
F%* is that we can find an IPM. More precisely 

Lemma 2.1. The process X e satisfies: 

1) For each function f € L X {D x Q,;P* D ) and t > 0: 

(13) E s *[f(X?,T X!/s u)]=M* D [fl 

2) For each function f € L l (dD x Q.;W% D ) and t > 0: 

(14) W*[( f(X*,T X s /s uj)dK E r ]=tM* dD [f}. 

Jo 

The first relation p3|) results from the structure of Cy (see 0), which has been defined so 
as to make e~ 2V ^dx invariant for the process X e . Once (|13p established, (jT^J) is derived from 
the fact that K e is the density of occupation time of the process X e at the boundary dD. 

2.3 Ergodic problems 

The next step is to determine the asymptotic behaviour as e — > of the quantities 

(15) / f{T X ,,^)dr and f /(r x . /e w) dK* r . 

Jo Jo 

The behaviour of each above quantity is related to the evolution of the process X e respectively 
inside the domain D and near the boundary dD. We shall see that both limits can be identified 
by solving ergodic problems associated to appropriate resolvent families. What concerns the first 
functional has already been investigated in the literature. The main novelty of the following 
section is the boundary ergodic problems associated to the second functional. 
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Ergodic problems associated to the diffusion process inside D 

First we have to figure out what happens when the process X e evolves inside the domain D. In 
that case, the pushing of the local time in ([T]) vanishes. The process X e is thus driven by the 
same effects as in the stationary case (without reflection term). The ergodic properties of the 
process inside D are therefore the same as in the classical situation. So we just sum up the main 
results and give references for further details. 

Notations: For p E [1; oo], L P (S1) denotes the standard space of p-th. power integrable functions 
(essentially bounded functions if p = oo) on (SI, Q, /i) and | • \ p the corresponding norm. If p = 2, 
the associated inner product is denoted by ( • , • ) 2 - The space C^°(D) (resp. C%°(D)) denotes 
the space of smooth functions on D with compact support in D (resp. D). □ 

Standard background: The operators on L 2 (Sl) defined by T x g(uj) — g(r x uj) form a strongly 
continuous group of unitary maps in L 2 (Sl). Let (ei,...,e,j) stand for the canonical basis of 
M. d . The group (T x ) x possesses d generators defined by Dig = limh^M^o h^ 1 (Th Gi g — g), for 
i = 1, . . . ,d, whenever the limit exists in the L 2 (il)-sense. The operators (Di)i are closed and 
densely defined. Given cp £ HiLi Dom(Di), Dip stands for the d-dimensional vector whose entries 
are Diip for i = 1, . . . , d. 

We point out that we distinguish Di from the usual differential operator d Xi acting on differ- 
entiable functions / : M. d — > M (more generally, for k > 2, d x .^ x . denotes the iterated operator 
d Xil ■ ■ ■ d Xik ). However, it is straightforward to check that, whenever a function ip E Dom(Z?) 
possesses differentiable trajectories (i.e. [i a.s. the mapping x h- > ip(r x oS) is differentiable in the 
classical sense), we have Diip(T x Lu) = d Xi ip(r x bj). 

We denote by C the dense subspace of 1? (fl) defined by 

(16) C = Span{g*ip;g E L°°(Q),<p £ C c °°(R d )} where g * <p{u) = [ g{ Tx uj)ip{x) dx 

Js. d 

Basically, C stands for the space of smooth functions on the random medium. We have C C 
Dom(Di) and Di(g -kip) = —g * d Xi tp for all 1 < i < d. This quantity is also equal to Dig ★ ip if 
g E Dom(A)- □ 

We associate to the operator C e (Eq. ([J}) an unbounded operator acting on C C L 2 (Vl) 

(17) L= -Di( aij D r ). 

Following [3 Ch. 3, Sect 3.] (see also [19l Sect. 4]), we can consider its Friedrich extension, still 
denoted by L, which is a self-adjoint operator on L 2 (Vl). The domain H of the corresponding 
Dirichlet form can be described as the closure of C with respect to the norm || vile = I VI2 + l-^Vli- 
Since L is self-adjoint, it also defines a resolvent family (U\)\ > o- For each / E L 2 (fl), the function 
w\ = U\(f) £ Ifl Dom(L) equivalently solves the L 2 (f2)-sense equation 

(18) \w x -Lw x = f 
or the weak formulation equation 

(19) VveH, X(wx,f)2 + (l/2)(a ij D i wx,D j <p) 2 = (f,<p) 2 . 
Moreover, the resolvent operator U\ satisifes the maximum principle: 

Lemma 2.2. For any function f £ L°°(n), the function U\{f) belongs to L°°{^1) and satisfies 

\Ux{f)U < l/loo/A. 

The ergodic properties of the operator L are summarized in the following proposition: 

Proposition 2.3. Given f E L (SI), the solution W\ of the resolvent equation Xw\ — Lw\ = f 
(X > ) satisfies 

\Xw x - M[/]| 2 ^0 as X -> 0, and VA > 0, |A 1/2 L>u; A |2 < A" 1/2 |/| 2 . 
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Boundary ergodic problems 

Second, we have to figure out what happens when the process hits the boundary dD. If we want 
to adapt the arguments in [22j . it seems natural to look at the unbounded operator in random 
medium _ff 7 , whose construction is formally the following: given w G £1 and a smooth function 
ip G C, let us denote by : D — > R the solution of the problem 

(2Q v f L^u^x) = 0, x G D, 

where the operator is defined by 

(21) L u f(x) = {l/2) aij {T x u)dl iXj f{x) + bi(T x u;)d x J(x) 
whenever / : D — > R is smooth enough, say / G C 2 (D). Then we define 

(22) H^(u)= li {u)d Xi u w {Q). 

Remark. Choose e = 1 in ([1]) and denote by (X ,K ) the solution of ([T]). The operator is 
actually the generator of the Q-valued Markov process Z t (uj) = Ty^^io, where Y t (uj) — X^^x^ 
and the function K^ 1 stands for the left inverse of K 1 : K^ 1 ^) — inf{s > 0; > t}. The 
process Z describes the environment as seen from the particle whenever the process X 1 hits the 
boundary dD. 



The main difficulty lies in constructing a unique solution of Problem (|20[) with suitable growth 
and integrability properties because of the lack of compactness of D. This point together with 
the lack of IPM are known as the major difficulties in homogenizing the Neumann problem in 
random media. We detail below the construction of H 1 through its resolvent family. In spite 
of its technical aspect, this contruction seems to be the right one because it exhibits a lack of 
stationarity along the ei-direction, which is intrinsec to the problem due to the pushing of the 
local time K e , and conserves the stationarity of the problem along all other directions. 

First we give a few notations before tackling the construction of H 1 . In what follows, the 
notation (xi,y) stands for a d-dimensional vector, where the first component x\ belongs to K 
(eventually K + = [0; +oo)) and the second component y belongs to R d_1 . To define an unbounded 
operator, we first need to determine the space that it acts on. As explained above, that space must 
exhibit a a lack of stationarity along the ei-direction and stationarity along all other directions. 
So the natural space to look at is the product space R+ x ft, denoted by equipped with the 

measure dfj, + "= dx\ ® d^i where dx\ is the Lebesgue measure on R + . We can then consider the 
standard spaces L p (fi + ) for p G [1; +oo]. 

Our strategy is to define the Dirichlet form associated to H 1 . To that purpose, we need to 
define a dense space of test functions on f2 + and a symmetric bilinear form acting on the test 
functions. It is natural to define the space of test functions by 

C(fi+) - Span{p(a!iMw);pG C~([0;+oo)), V G C}. 

Among the test functions we distinguish those that are vanishing on the boundary {0} x f2 of f2 + 

C c (r!+) = Span{p(a!i)v>(w); p G C c °°((0; +«>)), <p G C}. 

Before tackling the construction of the symmetric bilinear form, we also need to introduce 
some elementary tools of differential calculus on 57+ . For any g G C(fi + ), we introduce a sort 
of gradient dg of g. If g G C(£7 + ) takes on the form p(xi)(f(cu) for some p G C^°([0; +oo)) and 
ip G C, the entries of dg are given by 

d\g{xx,u) = d Xl g(xi)ip(u), and, for i = 2, . . . , d, dig(xi,uj) = p(x 1 )D i ip{uj). 
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We define on C(f2 + ) the norm 

(23) N(gf = \g(0,-)\l+ f \dg\ 2 3 d^+, 

Jn+ 

which is a sort of Sobolev norm on f2 + , and W 1 as the closure of C(£! + ) with respect to the norm 
N (W 1 is thus an analog of Sobolev spaces on f2 + ). Obviously, the mapping 

P : W 1 3 g h-> 0(0,-) £ L 2 (ft) 

is continuous (with norm equal to 1) and stands, in a way, for the trace operator on f2 + . Equip 
the topological dual space (W 1 )' of W 1 with the dual norm N'. The adjoint P* of P is given by 
P* : ip £ L 2 (f2) 1— » P*((f) £ (W 1 )' where the mapping P*(f) exactly matches 

P*( V ) : g e W 1 ^ (g, PV) = ( Vs fl(0, 0)2- 

To sum up, we have constructed a space of test functions C(f2 + ), which is dense in W 1 for the 
norm AT, and a trace operator on W 1 . 

We further stress that a function g £ W 1 satisfies dg — if and only if we have g(xi,u>) = f(u>) 
on fl + for some function / £ L 2 (Q) invariant under the translations {t x ;x £ {0} x For 
that reason, we introduce the cr-field Q* C Q generated by the subsets of fl that are invariant 
under the translations {t x ;x £ {0} x IR^ 1 }, and the conditional expectation Mi with respect to 
Q*. 

We now focus on the construction of the symmetric bilinear form and the resolvent family 
associated to H 1 . For each random function tp defined on fi, we associate a function ip + defined 
on f2 + by 

V(xi, u) £ f2 + , tp + (xi,u>) = ip(r Xl uj). 

Hence, we can associate to the random matrix a (defined in Section [T]) the corresponding matrix- 
valued function a + defined on f2 + . Then, for any A > 0, we define on W 1 x W 1 the following 
symmetric bilinear form 

(24) B x (g,h) = X(Pg,Ph) 2 + l- f at d i9 djh . 

1 Jn+ 

From ([5]), it is readily seen that it is continuous and coercive on W 1 x W 1 . From the Lax-Milgram 
theorem, it thus defines a continuous resolvent family G\ : (W 1 )' — > W 1 such that: 

(25) VFe (W 1 )', VgE W 1 , B x (G x F,g) = (g,F). 
For each A > 0, we then define the operator 

(26) R x ■■ L 2 {n) -> L 2 {n) 

if ^ PG x P*(ip) 
Given ip £ L 2 (il), we can plug F = P*cp into ([25]) and we get 

(27) VgeW 1 , B x (G x P*cp, g) = (g, P* v ), 
that is, by using (|24|) : 

(28) VgeW 1 , \{R xV ,Pg) 2 + l - I a± di{G x P*tp) d j9 dfx + = (g(0, ■), tp) 2 , 

2 Jn+ 

The following proposition summarizes the main properties of the operators (i?A)A>o> and in 
particular their ergodic properties: 
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Proposition 2.4. The family (R\)\ is a strongly continuous resolvent family, and: 

1) the operator R\ is self-adjoint. 

2) given ip G L 2 (fl) and A > 0, we have: 

if e Kei(\R x - I) cp = M 1 [ip]. 

3) for each function tp G L 2 (Vt), \\R\ip - Mi[tp]\ 2 — > as A — > 0. 

The remaining part of this section is concerned with the regularity properties of G\P*ip. 

Proposition 2.5. Given ip G C, i/ie trajectories of G\P*ip are smooth. More precisely, we can 
find N C f2 satisfying fi(N) = and smc/i f/ia£ Vw G O \ AT, i/ie function 

: x = (xi,y) G 5 h-> G\P*ip(xi, T( 0)J/ )w) 

belongs to C°°{D). Furthermore, it is a classical solution to the problem: 



(29) 



L u u u (x) = 0, x G D, 

Au w (a;) - 7 i (r x a;)S Zj {t u ,(x) = i/j^w), X G <9Z>. 



In particular, the above proposition proves that (R\)\ is the resolvent family associated to 
the operator H~ r This family also satisfies the maximum principle: 

Proposition 2.6. (Maximum principle). Given <p € C and A > 0, we have: 

|G A PVU~ (n +) < A- x bU. 



2.4 Ergodic theorems 

As already explained, the ergodic problems that we have solved in the previous section lead to 
establishing ergodic theorems for the process X e . The strategy of the proof is the following. First 
we work under P e * to use the existence of the IPM (see Section |2~2"|) . By adapting a classical 
scheme, we derive from Propositions 12.31 and 12.41 ergodic theorems under P £ * both for the process 
X e and for the local time K e : 

Theorem 2.7. For each function f G L 1 (fi) and T > 0, we have 



(30) 



lim E £ 

e^O 



sup 

0<*<T ' Jo 



f(r X e /e oj)dr-tU[f}\ 



= 0. 



Theorem 2.8. If f £ L 2 (fl), the following convergence holds 



(31) 



lim E £ 



sup 

0<t<T ' Jo 



/(rjt./ e w) tiff?- Mi | 



= 0. 



Finally we deduce that the above theorems remain valid under P e thanks to (jlljl . 

Theorem 2.9. 1) Let (f £ ) £ be a family converging towards f in L 1 (f2). For each fixed 6 > 
and T > 0, the following convergence holds 



(32) 



lim I 

e->0 



sup I / f e (Tx'/e")dr-tM[f] \ >5 

0<t<T Jo 



= 0. 



2) Let (f e ) s be a family converging towards f in L 2 (fl). For each fixed S > and T > 0, the 
following convergence holds 



(33) 



lim I 

e^O 



sup 

0<t<T ' Jo 



f e (T Xf/e u) dK e r - Mi[/]fff| > S 



0. 
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2.5 Construction of the correctors 



Even though we have established ergodic theorems, this is not enough to find the limit of equation 
([T]) because of the highly oscillating term e~ 1 b(r x ^ / £ oj) dt. To get rid of this term, the ideal 
situation is to find a stationary solution u 1 : — > R to the equation 

(34) - Lu l = b,. 

Then, by applying the Ito formula to the function u l , it is readily seen that the contribution of 
the term £~ 1 bi(r X z / E u>) dt formally reduces to a stochastic integral and a functional of the local 
time, the limits of which can be handled with the ergodic theorems 12.91 

The problem is that the lack of compactness of a random medium makes you cannot find a 
stationary solution to (f3~4")l . As already suggested in the literature, a good approach is to add 
some coercivity to the problem (|34p and define, for i = 1, . . . , d and A > 0, the solution u\ of the 
resolvent equation 

(35) Xu\ - Lu\ = bi. 

If we let A go to in of ([35]). the solution u\ should provide a good approximation of the solution 
of (p?4"]) . Actually, it is hopeless to prove the convergence of the family («\)a in some L p (f2)-space 
because, in great generality, there is no stationary L p (^-solution to (f3"l)) . However we can prove 
the convergence towards of the term Xu\ and the convergence of the gradients Du\: 

Proposition 2.10. There exists C € (L 2 (fl)) d such that 

(36) X\u\\l + \Du\ - C% -> 0, asX^O. 

As we shall see in Section [2^6l the above convergence is enough to carry out the homogenization 
procedure. The functions C (i < d) are involved in the expression of the coefficients of the 
homogenized equation ([6]). For that reason, we give some further qualitative description of these 
coefficients: 

Proposition 2.11. Define the random matrix- valued function £ G L 2 (0; R dxd ) by its entries 
dj — Cl — li m A^o DiU x . Define the matrix A and the d-dimensional vector T by 

(37) A = M[(I + C)a(l + C)], which also matches M[(I + C*)a], 

(38) f = M[(I + C*h] e M. d , 

where I denotes the d-dimensional identity matrix. Then A obeys the variational formula: 

(39) VA g R d , X*AX = inf M[(X + Dtp)* a(X + D(p)\. 

tp&C 

Moreover, we have A > AI (in the sense of symmetric nonnegative matrices ) and the first compo- 
nent Ti ofT satisfies T% > A. Finally, T coincides with the orthogonal projection Mi [(I + £*)-/]. 

In particular, we have established that the limiting equation ([6]) is not degenerate, namely 
that the diffusion coefficient A is invertible and that the pushing of the reflection term f along 
the normal to dD does not vanish. 

2.6 Homogenization 

Homogenizing ([1} consists in proving that the couple of processes (X e ,K e ) e converges as e — > 
(in the sense of Theorem II. 2 [) towards the couple of processes (X, K) solution of the RSDE |6]). 
We also remind the reader that, for the time being, we work with the function x( x ) = e~ 2V ( x \ 
We shall see thereafter how the general case follows. 

First we show that the family (A e ,AT e ) e is compact in some appropriate topological space. 
Let us introduce the space D([0, T}; R + ) of nonnegative right-continuous functions with left limits 
on [0,T] equipped with the S-topology of Jakubowski (see Appendix IF)) . The space C([0,T}; D) 
is equipped with the sup-norm topology. We have: 
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Proposition 2.12. Under the law P e , the family of processes (X £ ) £ is tight in C([0, T]; D), and 
the family of processes (K e ) £ is tight in D( [0, T]; K+). 

The main idea of the above result is originally due to Varadhan and is exposed in [161 Chap. 3] 
for stationary diffusions in random media. Roughly speaking, it combines exponential estimates 
for processes symmetric with respect to their IPM and the Garsia-Rodemich-Rumsey inequality. 
In our context, the pushing of the local time rises some further technical difficulties when the 
process X e evolves near the boundary. Briefly, our strategy to prove Proposition 12 . 1 21 consists in 
applying the method [161 Chap. 3] when the process X e evolves far from the boundary, say not 
closer to 3D than a fixed distance 9, to obtain a first class of tightness estimates. Obviously, these 
estimates depend on 9. That dependence takes place in a penalty term related to the constraint 
of evolving far from the boundary. Then we let 9 go to 0. The limit of the penalty term can 
be expressed in terms of the local time K e in such a way that we get tightness estimates for the 
whole process X e (wherever it evolves). Details are set out in the appendix [Gl 

It then remains to identify each possible weak limit of the family (X e , K e ) e . To that purpose, 
we introduce the corrector u x 6 L 2 (tt; R d ), the entries of which are given, for j = 1, . . . , d, by the 
solution u J x to the resolvent equation 

\u J x — Lu 3 x = bj . 

Let C £ L 2 (Q;'R dxd ) be defined by C,j = k m A^o DiU x (see Proposition 12. 10p . As explained in 
Section [^751 the function u x is used to get rid of the highly oscillating term e~ 1 b(r X e / £ co) dt in (TT|) 
by appling the Ito formula. Indeed, since /^-almost surely the function <j) : x i— > u x (t x lo) satisfies 
\<fr — L u <t> = b(r.Lo) on K d , the function x i— > u\(t x lo) is smooth (see [5] Th. 6.17]) and we can 
apply the Ito formula to the function x i— > eu\{T x / c uj). We obtain 

Edu x {T X ' t / e u) =-Lu x (t X i/ e u) dt + Du* x ~f(T X i/ E u) dK £ t + Du* x ct(tx- / e w) dB t 

(40) = \ <yXux ~ b ^ Tx t/^ dt + Du xl{ T Xt/e^) dK\ + Du* x cr{T X e/eU) dB t . 
By summing the relations (14T))) and ([1} and by setting A = e 2 , we deduce: 

(41) X\ =x - e(u e 2(T X c/ e uj) - u £ 2(t xs/£ u})) +e / u^(t x -- / e w) dr 

Jo 

+ [ (I + Du; 2 ) 1 (T X . /£ u)dK e r + [ (I + Du* 2 )a(T X s /£ cj)dB r . 
Jo Jo 

= x- G]' e + GY + Gf £ + M e t . 

So we make the term e b(r X e / e w) dt disappear at the price of modifying the stochastic integral 
and the integral with respect to the local time. By using Theorem 12.91 we should be able to 
identify their respective limits. The corrective terms G 1 ^ and G 2 - e should reduce to as e — * 0. 
This is the purpose of the following proposition: 

Proposition 2.13. For each subsequence of the family {X s , K £ ) e , we can extract a subsequence 
( still indexed with e > 0) such that: 

1 ) under P £ , the family of processes {X £ , M £ , K e ) £ converges in law in C([0, T];D)x C( [0, T] ; R d ) x 
-D([0, T]; M + ) towards (X, M, K), where M is a centered d-dimensional Brownian motion with co- 
variance 

A = M[(I + C*)a(I + C)] 

and K is a right- continuous increasing process. 

2) the finite- dimensional distributions of the families {G\' e ) £ , (G 2 ' £ ) e and (G 3,£ — TK e ) £ con- 
verge towards in W -probability, that is for each t 6 [0, T] 

V<5 > 0, limP £ (|G^| > 5) = (i = 1,2), lim P £ (|G?' £ - YKf\ ><j) = 0. 
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Proof. 1) The tightness of (X s , K £ ) results from Proposition 12. 121 To prove the tightness of the 
martingales (M £ ) e , it suffices to prove the tightness of the brackets (< M £ >) E , which are given 

by 

< M E > t = / (1 + Du* 2 )a(I + Du £ 2)(T X z/ £ uj)dr. 
Jo 

Proposition HHni and Theorem [231 lead to < M e > t -> At in probability in C([0, T];R dxd ) where 
A = M[(I + C)a(I + C)] • The martingales {M £ ) £ thus converge in law in C([0, T]; R d ) towards a 
centered Brownian motion with covariance matrix A (see [9]). 

2) Let us investigate the convergence of (G l,£ ) £ (i — 1,2). From the Cauchy-Schwarz inequality, 
Lemma I2TT1 and ([3S]) , we deduce: 



lim E £ 

e—0 



|£W^(^/ £ ^)| 2 + I / eu £ 2(T X , /£ co)dr\ 2 ] < (1 + i) hWe 2 ^ I2) = 0. 

We conclude with the help of ([TTjl . 

Finally we prove the convergence of (G 3 ' e ) e with the help of Theorem l2.9l Indeed, Proposition 
12.101 ensures the convergence of the family ((I + Du* 2 )j) e towards (I + C*)7 m L 2 (Vt) as e — > 0. 
Furthermore we know from Proposition ^, lll that T = Mi [(I+C*)"y] • The convergence follows. □ 

Since the convergence of each term in (|41j) is now established, it remains to identify the 
limiting equation. From Theorem IF. 21 we can find a countable subset S C [0,T[ such that the 
finite-dimensional distributions of the process (X £ , M £ , K e ) £ converge along [0, T] \ S. So we can 
pass to the limit in (HTj) along s, t € [0, T] \ S (s < t), and this leads to 

(42) X t =X s +A l ' 2 {B t -B s )+f{K t -K s ). 

Since (|42|) is valid for s,t £ [0,T] \ S (note that this set is dense and contains T) and since 
the processes are at least right continuous, (|l2")) remains valid on the whole interval [0, T]. As a 
by-product, K is continuous and the convergence of (X £ , M £ , K £ ) £ actually holds in the space 
C([0,T];L>) x C([0,T];R d ) x C([0,T];R + ) (see LemmaES}- 

It remains to prove that K is associated to X in the sense of the Skorokhod problem, that is 
to establish that {Points of increase of K} C {t; X\ = 0} or / Q X Y r dK r = 0. This results from 

the fact that Ve > X^ £ dK £ = and Lemma fF.41 Since uniqueness in law holds for the 
solution (X, K) of Equation (|42|) (see [23]), we have proved that each converging subsequence 
of the family (X £ ,K £ ) £ converges in law in C([0,T]; D x R + ) as e — > towards the same limit 
(the unique solution (X,K) of ©). As a consequence, under P e , the whole sequence (A £ ,X e ) e 
converges in law towards the couple (X, K) solution of (HJ). 



Replication method 

Let us use the shorthands Cd and C+ to denote the spaces C([0, T], L>) and C([0, T],R+) respec- 
tively. Let E denote the expectation with respect to the law P of the process (X, K) solving the 
RSDE © with initial distribution P(X € dx) = e- 2V{x) dx. From [23], the law P coincides with 
the averaged law f' x {-)e~ 2V ^dx where P^ denotes the law of (X, K) solving (|42[) and starting 
from x G D. 

We sum up the results obtained previously. We have proved the convergence, as e — > 0, of 
E £ [F(X £ ,K £ )\ towards E[F(X,K)}, for each continuous bounded function F : C D x C+ -> R. 
This convergence result is often called annealed because E e is the averaging of the law P£ with 
respect to the probability measure P* D . 

In the classical framework of Brownian motion driven SDE in random media (i.e. without 
reflection term in Q), it is plain to see that the annealed convergence of X s towards a Brownian 
motion implies that, in PJ,-probability, the law P| of X s converges towards that of a Brownian 
motion. To put it simply, we can drop the averaging with respect to ¥* D to obtain a convergence in 
probability, which is a stronger result. Indeed, the convergence in law towards of the correctors 
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(by analogy, the terms G 1:£ ,G 2,£ in (|4ip ) implies their convergence in probability towards 0. 
Moreover the convergence in P^-probability of the law of the martingale term M e in (|4"Tj) is 
obvious since we can apply [9] for P^-almost every (x,ui) G D x fi. 

In our case, the additional term G 3,e puts an end to that simplicity: this term converges, 
under the annealed law P £ , towards a random variable TK, but there is no obvious way to switch 
annealed convergence for convergence in probability. That is the purpose of the computations 
below. 

Remark and open problem. The above remark also raises the open problem of proving a so- 
called quenched homogenization result, that is to prove the convergence of X e towards a reflected 
Brownian motion for almost every realization u> of the environment and every starting point 
x G D. The same arguments as above show that a quenched result should be much more difficult 
than in the stationary case 



So we have to establish the convergence in P^-probability of E X [F(X £ , K £ )] towards E X [F(X, K)] 
for each continuous bounded function F : Cd x G + — > R. Obviously, it is enough to prove the 
convergence of E%[F(X £ , K £ )] towards E X [F{X, K)\ in L 2 (D x f2, P* D ). By using a specific feature 
of Hilbert spaces, the convergence is established if we can prove the convergence of the norms 



(43) M* D (E X [F(X £ ,K £ )]Y — > M* D (E X [F(X,K)}) 



as e — * 0, 



as well as the weak convergence. Actually we only need to establish (|4"3"|) because the weak 
convergence results from Section [2T6l as soon as (|43| is established. 

The following method is called replication technique because the above quadratic mean can 
be thought as of the mean of two independent copies of the couple (X e ,K s ). We consider 2 
independent Brownian motions {B l ,B 2 ) and solve (fT]) for each Brownian motion. This provides 
two independant (with respect to the randomness of the Brownian motion) couples of processes 
(X s ' 1 , K 6 ' 1 ) and (X £ > 2 , K e > 2 ). Furthermore, we have 



W x [F(X*, K £ )]f\ = U* D [E xx [F(X s >\K s > l )F(X^,K*> 2 )]] 

where E xx denotes the expectation with respect to the law P^ of the process (X 6 ' 1 , K £:1 ,X 5,2 , K e ' 2 ) 
when both X s ' 1 and X e ' 2 start from x G D. Under M* D ¥ e xx , the results of subsections [231 [2751 
and Proposition 12.121 remain valid since the marginal laws of each couple of processes coin- 
cide with ¥ £ x . So we can repeat the arguments of subsection 12.61 and prove that the processes 
(X 6,1 , K 6 - 1 , X £ ' 2 , K e - 2 ) e converge in law in Cd x C + x Cd X D + , under M* D E% X , towards a process 
(X 1 , K 1 , X 2 ,K 2 ) satisfying: 



(44) Mt e [0, T], X\ = Xl + A 1/2 B] + YK], X 2 = X 2 + A 1/2 B 2 + TK, 



1 1 



where (B X ,B 2 ) is a standard 2d-dimensional Brownian motion and K X ,K 2 are the local times 
respectively associated to X 1 ,^ 2 . Let P denote the law of (X 1 , K 1 , X 2 , K 2 ) with initial distri- 
bution given by P(X X e dx,X 2 G dy)_= 5 x {dy)e~ 2V dx _and f xx the law of (X 1 , K 1 , X 2 , K 2 ) 
solution of pi)) where both X 1 and X 2 start from x G D. To obtain (|4"3"|) , it just remains to 
remark that 

I[F(X 1 ,K 1 )F(A >2 ,i? 2 )] = f E XX [F{X\ K^FiX 2 , K 2 )]e~ 2V ^dx 

Jd 

E x [F(X\ K 1 )]^ [F(X 2 , K 2 )] e~ 2V ^dx, 



D 



since, under ¥ xx , the couples (X 1 ^ 1 ) and (X 2 ,K 2 ) are adapted to the nitrations generated 
respectively by B 1 and B 2 and are therefore independent. □ 
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2.7 Conclusion 



We have proved Theorem 1 1.21 for any function x that can be rewritten as x( x ) — e~ 2V ( x \ where 
V : D — > R is defined in (fl"2")) . It is then plain to see that Theorem 11.21 holds for any nonnegative 
function \ not greater than Ce~ 2V ( x \ for some positive constant C and some function V of the 
type (112[) . Theorem 11.21 thus holds for any continuous function \ with compact support over D. 

Consider now a generic function x '■ D — > R+ satisfying Jg x( a; ) das = 1 and x' : D — > R + with 
compact support in £>. Let j4 £ C x i) be defined as 

A e = {(u,x) efixfl; |E|(F(X £ ,i!r E )) -E X (F(X,#))| > 6} . 

From the relation M J D l A sx(x)dx <MJ D \x(x) - x'(x)\dx +MJ D l A ex'(x)dx, we deduce 

limsupM / l A ex{x)dx <M \x(x) - x'(x)\dx, 

e-»0 JD J D 

in such a way that the Theorem 11.21 holds for x by density arguments. The proof is completed. 

Appendix 

A Preliminary results 

Notations: Classical spaces. Given an open domain O C R n and k e N U {oo}, C k (0) (resp. 
C k ((D), resp. C k {0)) denotes the space of functions admitting continuous derivatives up to 
order k over O (resp. over O, resp. with continuous bounded derivatives over D). The spaces 
C k (0) and C k {0) denote the subspaces of C k (0) whose functions respectively have a compact 
support in O or have a compact support in O. Let C\' 2 denote the space of bounded functions 
/ : [0, T] x D — > R admitting bounded and continuous derivatives dtf, d x f, d 2 x f and d xx f on 
[0, T] x D. 

Green's formula: 

We remind the reader of the Green formula (see [T5J eq. 6.5]). We consider the following operator 
acting on C 2 (D) 

p 2V(x) d 

(45) £ s v = — ^— 9 Xi (e- 2V ^ aij (r x/e u)d Xj ) , 

where V : D — > R is smooth. For any couple (ip, ip) £ C 2 (D) x C\{D), we have 

/ C E v <p(x)ip(x)e- 2V Wdx + ± [ a l3 (T x/£ Lu)d x Mx)d x M x ) e ~ 2V(x) dx 

J D * JD 

(46) =-\( li {r x/e u)d x M^{x)e- 2V ^ dx. 

1 JdD 

Note that the Lebesgue measure on D or 3D is indistinctly denoted by dx since the domain of 
integration avoids confusion. 

PDE results: 

We also state some preliminary PDE results that we shall need in the forthcoming proofs: 

Lemma A.l. For any functions / € C^°(D) and g,h G C£°(D), there exists a unique classical 

solution w £ G C°°([0, T];D) D C^' 2 to the problem 

(47) 

d t w e — £yw e + gw e + h on [0, T] x D, •y i (r./ e u)d Xi w e = on [0, T] x 3D, and w e (0, •) = /. 
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Proof. First of all, we remind the reader that all the coefficients involved in the operator Cy 
belong to C£°(D). From [13l Th V.7.4], we can find a unique generalized solution w' e in Cj/ 2 to 
the equation 

d t w' E = C v w' E + gw' £ + C v f + gf + h, w' e (0, •) = on D, j(T. /e uj)d Xt w' £ = on [0, T] x dD. 
From [13 1 IV. 10], we can prove that w' e is smooth up to the boundary. Then the function 

w E (t, x) = w' £ (t, x) + f(x) £ C°°([0, T] x D) n Cl' 2 
is a classical solution to the problem (j47f . □ 



Lemma A. 2. The solution w £ given by Lemma \A.l\ admits the following probabilistic represen- 
tation: V(i, x) £ [0, T] x D, 



t , ft 



w £ (t,x)=W x *[f(X e t )vq>(j o g{X £ )dr) + J h{X £ ) exp ( J g(X s u )du)dr 

Proof. The proof relies on the Ito formula (see for instance [TQl Ch. II, Th. 5.1] or 6, Ch. 2, 
Th. 5.1]). It must be applied to the function (r, x, y) i— > w e (t — r, x)exp(y) and to the triple of 
processes (r, X £ , f£ g(X £ )du). Since it is a quite classical exercise, we let the reader check the 
details. □ 

B Proofs of subsection 12.21 

Proof of Lemma \2.1\ 1) Fix t > 0. First we suppose that we are given a deterministic function 
/ : D — > R belonging to C%°(D). From Lemma [A. 11 there exists a classical bounded solution 
w £ e C°°([0, t]xD)n C^ 2 to the problem 

d t w e = CyW £ on [0,i] x D, 7 i (r./ e w)9 a;i ?« e = on [0,t] x dD, and w e (0, ■) = /(■), 

where Cy is defined in (|45[) . Moreover, Lemma lA. 21 provides the probabilistic representation: 

«; e (t,a:)=E|*[/(X t e )]. 

The Green formula ([4^)) then yields 

ft / u) £ (t,x)e- 2V(2;) da ; = / C £ v w e {t,x)e- 2V{x) dx 
Jd Jd 

= -\ [^ t (r x /^)d Xt w e (t,x)e- 2V ^ dx = 



dD 



so that 



(48) / W*[f{X e )]e~ 2V{x) dx= / f{x)e- 2V(x) dx. 

Jd Jd 



It is readily seen that (|48[) also holds if we only assume that / is a bounded and continuous 
function over D: it suffices to consider a sequence (f n )n C C^ D {D) converging pointwise towards 
/ over D. Since / is bounded, we can assume that the sequence is uniformly bounded with respect 
to the sup- norm over D. Since (|48p holds for f n , it just remain to pass to the limit as n — > oo 
and apply the Lebesgue dominated convergence theorem. 

We have proved that the measure e~ 2VtyX ^ dx is invariant for the Markov process X s (under 
P e *). Its semi-group thus uniquely extends to a contraction semi-group on L 1 (L>, e~ 2V ^ dx). 
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Consider now / £ L X (D x £l;¥* D ) and e > 0. Then, \i almost surely, the mapping x \— * 
f(x,T x / e u>) belongs to L}(D, e~ 2V ^ dx). Applying (|35)) yields, fi almost surely, 

\l*[f{Xlr xtJe u)]e- 2V ^dx= [ f{x,T x/e u)e- 2V ^dx. 

D J D 

It just remains to integrate with respect to the measure fi and use the invariance of /j, under 
translations. 

Let us now focus on the second assertion. As previously, it suffices to establish 
E r [ / fi X r) dK r] e- 2V{x) dx = t [ f(x)e- 2V ^ dx 

D JO JdD 

for some bounded continuous function / : dD — > R. We can find a bounded continuous function 
/ : D — ► R such that the restriction to dD coincides with / (choose for instance / = / o p where 
p : D — > dD is the orthogonal projection along the first axis of coordinates). 

Recall now that the local time Kf is the density of occupation time at dD (see [H Prop. 1.19] 
with tp(x) = X\, Vo = 7 and a 2 (x) = 1). Hence, by using ([35]). 

¥.%*[[ f(X e r )dK e r ]e~ 2V{x) dx= ( W* [limT 1 f f(X E r )l m (X^ E ) dr] e~ 2V W dx 



D JO JD "5->0 j q 



= lim / Ef [S- 1 / f(Xt)l m (X^)dr\e- 2V W dx 
°^°Jd Jo 

^thmS- 1 { f(x)l [QS ](xi)e~ 2V{x) dx 
=t f f(x)e~ 2V(x) dx. □ 

JdD 

C Proofs of subsection 12.31 

Generator on the random medium associated to the diffusion process inside D 

Proof of Provosition WJA The first statement is a particular case, for instance, of [191 Lemma 6.2]. 
To follow the proof in [19], omit the dependency on the parameter y, take H = and = f. 
To prove the second statement, choose cp = w\ in (TT9"]) and plug the relation 

(f,w x ) 2 < \f\ 2 \w x \ 2 < 1/(2A)|/H + (A/2)|«; A || 

into the right-hand side to obtain AI^aII + { a ijDiW\, DjW\J < |/||/A. From (JSJ), we deduce 
A\Dw x \l < l/ll/A and the result follows. □ 



Proof of Lemma \2.2\ The proof is quite similar to that of Proposition 12.61 below. So we let the 
reader check the details. □ 

Generator on the random medium associated to the reflection term 

Proof of Proposition (|2.4| . The resolvent properties of the family (R\)\ are readily derived from 
those of the family (G\)\. 

So we first prove 1). Consider tp, tp € L 2 (n). Then, by using (|26|) and (I27j) . we obtain 

(Rx<p,iPh =( p GxP*f^)2 = (G x P*<p,P*iP) = B x (G x P*4>,G x P*ip) 
=B x (GxP* i p 1 GxP*rp) = (G A P>,P» = (v,^a^) 2 

so that R\ is self-adjoint in L 2 (il). 
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We now prove 2). Consider <p £ L 2 (f2) satisfying \R x <p = <p for some A > 0. We plug 
= G x P*ip e W 1 into 



(49) \\Rxv\l + \J n+ a±d i (G x P*<p)d j (G x P*<p)dn + = (PG x P*<p, V ) = (R x <p,<p) 2 . 

Since XR x cp = tp, the right-hand side matches (R x tp,tp)2 = A|i?A¥>|2 so that the integral term 
in (|49|) must vanish, that is / n+ di(G x P*tp) dj(G\P*tp) dfi + = 0. From j5]), we deduce 
d(G x P*tp) = 0. Thus, GaFV(0,-) is 5* -measurable. Moreover, we have \G x P*ip(0,-) = 
XPG x P*tp = \R x tp = tp so that tp is Q* -measurable. Hence tp = Mi[y>]. 

Conversely, we assume tp — Mi[cp], which equivalently means that tp is Q* measurable. We 
define the function u : £1+ — ► M by u{x\,uS) = <p(u). It is obvious to check that it belongs to W 1 
and satisfies du — 0. So B x (u, •) = (•, XP*(p) for any A > 0. This means it = XG x P*(p in such a 
way that XR x tp = \PG x P*tp = P{\G x P*tp) = Pu = tp. 

We prove 3). Consider <p € L 2 (il). Since the relation (|49l) is valid in great generality, (|49|) 
remains valid for such a function <p. Since the integral term in (|4"9")1 is nonnegative, we deduce 
A|-RavII - (R\ l P, l P)2 < \R\(p\2\<p\2- Hence |Ai?AV|2 < Ivh for any A > 0. So the family 
(XR x tp) x is bounded in L 2 (£l) and we can extract a subsequence, still indexed by A > 0, such 
that (XR x tp) x weakly converges in L 2 (Q) towards a function tp. Our purpose is now to establish 
that there is a unique possible weak limit tp = Mi[y] for the family (\R x tp) x . 

By multiplying the resolvent relation (A — /j,)R x R^tp = R^tp — R x tp by /i and passing to 
the limit as /i — > 0, we get \R x tp = (p. This latter relation implies (see above) that tp is im- 
measurable. To prove tp — Mi [</?], it just remains to establish the relation (tp, 1/1)2 = (tp, 1/1)2 
for every t/* -measurable function ip e L 2 (J7). So we consider such a function Obviously, it 
satisfies the relations Mj [(/)] = iji and \R x i/i = 1J1 (see the above item 2). We deduce 

{<p, 1/1)2 = {<p,\R\i/>)2 = \im(\R x ip,i{>)2 = {tp, 1/1)2- 

A — >0 

As a consequence, we have tp = Mi[y] and there is a unique possible limit for each weakly 
converging subsequence of the family (XR x tp) x . The whole family is therefore weakly converging 
in L 2 {&). 

To establish the strong convergence, it suffices to prove the convergence of the norms. As a 
weak limit, ip satisfies the property \<p\2 < liminfA^o \^R\<p\2- Conversely, (|4*9")) yields 

limsup|Ai?A¥'|2 < hnisup(Ai?A¥>, ¥>)2 = {&,<?) 2 = \<p\l 

A^O A^O 

and the strong convergence follows. □ 
The remaining part of this section is concerned with the regularity properties of the operator 
G X P* (Propositions [23] and [2?6]) and may be omitted upon the first reading. Indeed, though they 
may appear a bit tedious, they are a direct adatation of existing results for the corresponding 
operators defined on D (not on f2 + ). However, since we cannot quote proper references, we give 
the details. 

Given u 6 L 2 (f2 + ), we shall say that u is a weakly differentiable if, for i = 1, . . .,d, we can 
find some function diU £ L 2 (f2 + ) such that, for any g £ C c (f2 + ): 

udigd[i + = — I diugdfj, + . 

It is straightforward to check that a function u £ W 1 is weakly differentiable. For k > 2, the 
space W fe is recursively defined as the set of functions u S W 1 such that <9;it is fc— 1 times weakly 
differentiable for i = 1, . . . , d. 



Proposition C.l. If cp belongs to C, then G x P*tp £ flfeLi w 



A- 
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Proof of Provosition \C.1[ The strategy is based on the well-known method of difference quotients. 
Our proof, adapted to the context of random media, is based on Sect. 7.11 & Th. 8.8]. The 
properties of difference quotients in random media are summarized below (see e.g. |19[ Sect. 5]): 

i) for j = 2, . . . , d, r £ M\ {0} and g £ C c (f2+), we define 

A J r g(xi,uj) = -{g(xx,T rej bj) - g(xi,w)). 

ii) for each r£l \ {0} and g £ C c (f2+), we define 

Ka= ~(g(xx + r,uj) - g(xi,u)). 
r 

iii) for any j = 1, . . . , d, r € K\ {0} and g,h £ C c (0 + ), the discrete integration by parts holds 

/ A^gh dfi + = - gA j _ r hd/j, + 
Jn+ Jn+ 

provided that r is small enough to ensure that A J r g and A J r h belong to <D c (r2+). 

iv) for any j = 1, . . . , d, r £ M \ {0} and g £ <D c (f2+) such that A J r g £ C c (fi+), we have 



/ |A^| 2 d/i+< / \d j9 \ 



d^ 



Up to the end of the proof, the function G\P*cp is denoted by u. The strategy consists 
in differentiating the resolvent equation B\(u,-) = (-,P*tp) to prove that the derivatives of u 
equations of the same type. For p = 2, . . . , d, it raises no difficulty to adapt the method explained 
in [HO Sect. 5] and prove that the "tangential derivatives" d p u belongs to W 1 and solves the 
equation 

(50) B x (d p u,-) = (-,P*D pV >)-F p (-), 

where F p : W 1 — > R is defined by 



F p (g) = (1/2) f D p a± d t ud 3 gd^ 



In particular, dijU £ L 2 (Q + ; fi + ) for (i, j) ^ (1, 1). We let the reader check the details. 

The main difficulty lies in the "normal derivative" d\u: we have to prove that d\U is weakly 
differentiable. Actually, it just remains to prove that there exists a function d^u £ L 2 (f2 + ;/x + ) 
such that \/g £ C c (f2+): 



(51) / df 1 ugdfi + = — / d\udigdn + . 

To that purpose, we plug a generic function g £ C c (f2+) into the resolvent equation ([28]) . The 
boundary terms (P*ip,g) = (<p,g(0, -))2 and X(Pu, Pg)2 vanish and we obtain: 



d „ 

^ / a+ j d i ud j gdfi + = 0. 



We isolate the term corresponding to i = 1 and j = 1 to obtain (remind that an = 1) 



/ diudigdfi + =- ^2 / a+ j d i ud j gdfi + 

J + d j a±d i ugdfj, + + ^ J + a tj d ij u 9d^ 
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Since 8yu € L 2 (fi + ;^i + ) for (i, j) 7^ (1, 1), we deduce that 



/ diudigdfi + <C( g 
Jn+ JQ+ 



2 ^ + ) 1/2 



for some positive constant C. So the mapping g G C c (f2+) 1— ► J n+ d\udig dfi + is L 2 (fl + ; /i + )- 
continuous and there exists a unique function denoted by d\ x u such that (f5Tj) holds. As a conse- 
quence, <9iit is weakly differentiable, that is u G W 2 . Note that ([50)1 only involves the functions 
a, y and their derivatives in such a way that we can iterate the argument in differentiating (|50[) 
and so on. So it is clear that the proof can be completed recursively. □ 

Proof of Provosition [K5l The function u still stands for G\P*ip. From Proposition lC.li we have 
u G flfcli W fe an d it is plain to deduce that /i a.s. the trajectories of u are smooth and 

(52) Vx=(x 1) y)eD, d Xi u u (x) = d^x^T^u)). 

We let the reader check that point (it is a straightforward adaptation of the fact that an infinitely 
weakly differentiable function / : D — * R is smooth). 

It remains to prove that solves (|29p . To begin with, we state the following lemma 

Lemma C.2. For each function v € W 1 , we define : (xi,y) G D t— > v(x\,T^ y )Lo) . Then for 
every g G C^°(D) and ip G C we have: 



M[*P(lj) / (Xv u (y)-j i (T y Lj)d Xi v u (y))g{y)dy] = B x (v,ip* g) + M[if>(w) / L u v w (x)g{x) dx] 
JdD 

where the function ip * g : il + — > R belongs to W 1 and is defined by: 

tp*g(xi,u)= g(xi,-y)ij)(TyUj)dy. 
Let us consider g G C^°(D), ip G C. We first point out that 



B x (u, xJj * g) = (i/> * g, P*<p) = M[(p(u) / xP(T {0 . y) uj)g(0, -y) dy] = M[i/>(u>) / (p(T y u)g(y) dy] . 

jRrf-l JdD 

Then, by using Lemma lC.2l and the above relation, we obtain 

M[tp{u) / (Au w (y)-7i(Tj/w)9 x< {i w (2/)-v5(r^a;))e(y)rfy] =M[^(w) / L u u u (x)g(x) dx 

JdD Jd 

Since the above relation is valid for any ip G C, we deduce that /1 a.s. we have 

(Am w (2/) ~ fi(T y uj)d Xi u L0 (y) - tp{T y u))g{y) dy = L^u^(x)g(x) dx. 
dD Jd 

By choosing in turn a generic function g vanishing or not on the boundary, we deduce that fi a.s. 
we have: L^u^ = on D and Xu UJ (y) — j i (T y aj)d Xi u UJ (y) = ip(T y ui) for y G dD. □ 

Proof of Lemma \C.S\ First apply the Green formula (|46|) (with V = and e = I ): 

(Auo, (2/) -Ji(r v u>)d Xi v u ,(y))g(y)dy = / Xv u (y)g(y) dy + - f a ij (T x u>)d Xi v u ,(x)d X:j p(x)dx 

dD JdD 1 Jd 

L u v u (x)g(x) dx. 

Then we multiply the above relation by tp and integrate with respect to M. By using the invariance 
of \x under translations, we have 

M[VH / Xv u (y)g(y)dy] =XM[ / ^(uj)v{0,T y iu)g(y) dy] = XM[v(0,u) [ iP(T_ v uj)g{y) dy] 

JdD JdD JdD 

=X{Pv,Pip* g) 2 . 
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With similar arguments and (|52p . we prove 

i / a i j(T x u))dx i v u (x)d Xj p(x)dx =M / a+divdjip * gdp + . 

2 JD JR + 

The lemma follows. □ 



Proof of Proposition \2.6l We adapt the Stampacchia truncation method. More precisely, we 
introduce a function H : R — > R of class C 1 (R) such that 

i)Vs € R, |-ff'(s)| < C, m)Vs > 0, H'(s) > 0, m)Vs < 0, H'(s) = 0. 

We define K = \<p\oa/X and u x = G\P*tp. We let the reader check that H(u x - K) € W l . Then 
we plug g = H(u\ — K) into (|2"5|) and we obtain: 

A(Pm a , PiJ(u A - K)) 2 + \ [ a+ d t u x d ]Ux H'(u x - K) dp + = (PH(u x - K), ip) 2 . 

By subtracting the term \(K, H{Pu x — K)) 2 in each side of the above equality, we obtain: 

\{Pu x ~ K, H(Pu x - K)) 2 + I f a± 8 lUx d 3 u x H'{u x - K) dp+ = (H(Pu x -K),<p- \K) 2 . 

2 Jn+ 

Observe that the right-hand side is negative since cp — XK < and H(s) > for any s > 0. 
Furthermore, the left-hand side is positive since H'(s) > and sH(s) > for s E R. We deduce 
that both terms of the left-hand side reduce to 0. The relation X(Pu x — K, H(Pu x — K)) 2 = 
and the properties of H (sH(s) > for sel and sH(s) > for s > 0) ensure that Pu x — K < 0, 
that is 

(53) PH(u x - K) = 0. 

The relation \ J n+ diU X djU X H'(u x — K) dfj, + — and © prove that |9ua| 2 H'(u x - K) = 
p + a.s.. In particular, we deduce that 

(54) d(H(u x - K)) = 0. 

By gathering {53]) and we deduce N(H(u x - K)) = (recall the definition of N in (|2"3"|) ). 
So H(u x — K) = and this means u x < K. □ 

D Proofs of subsection 12.41 

Proof of Theorem \2. 7[ We first suppose that / belongs to C. Even if it means replacing / by 
/ — M[/], it is enough to treat the case M[/] = 0. We consider the solution v x g L 2 (f2)nDom(£) 
to the resolvent equation 

(55) Xv x - Lv x = f. 

For the same reason as in the proof of Proposition ^. 51 p a.s. the function i3 : x £ R d i— > v x (r x u>) 
satisfies Ai?(x)-i w ^(a;) = f(r x oj) x e R d . So -d is smooth Th. 6.17]. Applying the Ito formula 
to the function x <— > v x (t x lo) then yields 

dv\(T X */ e u) =£~ X DiV\(Tij(Txt /e^dBf 3 - E~ x d Xi V(X%)a,ij DjVx(rxf/ e u) dt 
+ E- 2 Lv x (T X z/ e u) dt + s~ 1 D l v x ~f t (Tx-/ e u) dK e t . 
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In the above expression, we replace Lv x by A^a — f, multiply both sides of the equality by s 2 
and isolate the term f(T X e/ e L>) dt. We obtain 

\ f(r X f/e^)dr=£ D l v X (T lj {T Xc j e uj)dB* r : i ~e 2 {vx{Tx' i /e^)-v x {T X! j e u:))+ i \v x {T Xf/e uj)dr 
Jo Jo Jo 



(56) 



e DtVxJ^Txf/e^dK^-e d Xl V(X^.)a lj D j v x (T X , /£ uj)dr 
Jo Jo 



def A l,s,A _ A 2,e,A + A 3,e,A + A 4,e,A _ A 5,e,A^ 



Let us investigate the quantities A 1,£ > A , A 2 < e ' A , A 3 ' e ' A , A 4 ' £ ' A and A 5 ' £ ' A . Using the Doob in- 
equality and Lemma |2. 11 we have: 

E e *[ sup \Al^ x \ 2 } <4e 2 TM* D [\D z v x( T t3 \ 2 } <Ce 2 \Dv x \ 2 2 

0<t<T 

for some positive constant C only depending on T and |<r|oo- Hence E e * [sup < t <r |Aj' e,A | 2 ] — > 
as e — > 0, for each fixed A > 0. Similarly, by using the boundedness of a, 7, d x V, we can prove 

E £ *[ sup |A^ e ' A | + sup |A^ e ' A | 2 ] ^0, ase^O. 

0<t<T 0<t<T 

From Lemma v x is bounded by |/|oo/A. We deduce 

W*[ sup |A 2 < £ < A | 2 ] < 4£ 4 |/|^A- 2 -> 0, ase^O. 

0<t<T 

By taking the limsup e _, in (|56p and by using the convergences of A 1,e ' A , A 2 ' e ' A , A 4,e ' A , A 5,E ' A 
towards 0, we deduce 

limsupE E *[ sup I / f(T X e /e u)dr\] < limsupE e *[ sup \Af e ' x \]. 

e-*0 0<t<T Jo T e^O 0<t<T 

Furthermore, from Lemma |2. 11 we have 

limsupE"[ sup \A 3 t ' £ ' X \] < limsup / E" [\Xv x (t xc /£ uj)\] dr = T\Xv x \i < T\Xv 

e^O 0<t<T e^O Jo 



A 2- 



From Proposition 12.31 we have IA^a^ — > as A goes to 0. So it just remains to choose A small 
enough to complete the proof in the case of a smooth function f E C. The general case follows 
from the density of C in and Lemma [2~T1 □ 



Proof of Theorem \2.8[ Once again, from Lemma 12.11 and density arguments, it is sufficient to 
consider the case of a smooth function f E C. Even if it means replacing / with / — Mi[/], 
it is enough to consider the case Mi[/] = 0. Let us define, for any A > 0, u x = G x P*f and 
fx = R\fi the definitions of which are given in Section l2~3l (boundary ergodic problems). We still 
use the notation u x (x) = u x (xi, t^q^oj) for any x — (xi, y) € D. We remind the reader that the 
main regularity properties of the function {t A are summarized in Proposition 12.51 In particular, 
/i a.s., the mapping x i— > u x (x) is smooth and we can apply the ltd formula: 

d{eut{XI/e)) = [e^L^utiX! /e) - d Xj V(X!) aiJ (T X . /£ Lo)d Xz ut(X! /e)} dt 
(57) + d^uliXt/e)^^;^) dB? + ll {T XU£ oj)d Xi u x jXt/e) dK\ 

In the above expression, we use the relation L^u^ = inside D. Furthermore, since 1id Xi u x (x) = 
^f x { T xu) — f{r x uj) on dD and dK\ = l.gr)(Xf)dK%, we deduce 

l i {r xt/ ^)d Xi u x U3 {XI/e)dKt = (Xf x - f)(r Xf/e cj) dK\. 



22 



Hence, ((57)) yields 

f(r X ; /e w) dK = - (eut(X;/e) - £u x u {X s Jej) - f d X} V(X e r ) aij (T X s /e u;)d Xi ut(X e r /e) dr 



Jo Jo 

(58) = - A]' 5 A 2 ' £ + A?' £ + Af s . 

The next step of the proof is to prove that A ,e , A 2,e , A 3 ' e converge to as £ goes to for each 
fixed A > 0. Clearly, from Proposition 12.61 we have 

W*[ sup |Aj' B | 2 ] < 4e 2 \u x \ 2 Laa{n+) — >0. 

Let us now focus on A 2 ' £ . We use the boundedness of d Xj V, (1 <i,j< d) and Lemma |2~H 
W*[ sup |A 2 > E | 2 ] <T\a\lx S up\d x V\ 2 xM* D [\d x ut(-/e)\ 2 ]. 

0<t<T D 

Furthermore 

M* D [\d x u x J : /e)\ 2 ]=m [ \du x (x 1 /e,r y/e w)\ 2 e- 2V ( x ^dx 1 dy 

(59) =m[ \du x { Xl /e,uj)\ 2 ( [ e ~ 2V(xi < v) dy)d Xl . 

We point out that the function V given by (fT2")l satisfies 

d - f f „-2V(xt,v) , 



(60) S = sup / e- 2V ^> dy <+oo. 

By gathering (|60[) and (f59"|) and by making the change of variables u = xi/e, we deduce that 
M* D [\d x u*(-/s)\ 2 ] is not greater than eSj n+ \du x \ 2 dfi+. So W* [sup < t < T |A 2 ' e | 2 ] converges 
to as e — ► 0. By combining the same argument with the Doob inequality, we prove that 

E"[sup < f < T |A^| 2 ] ->0ase->0. 

So, taking the limsup e _ >0 in (|58p and using the above convergences yields 

limsupE e *[ sup | f f{T X e /e w)dK s r \] <limsu P E £ *[ f \Xf x {t x , /e u)\ dK$] . 

e^O 0<t<T Jo e^O Jo 

By using Lemma l2.1l in the right-hand side of the previous inequality, we deduce, for any A > 0, 
limsupE"[ sup | [ f(T X e /e uj)dK e r \] < TM* dD [\Xf x \] = \\f x \ t T f e~ 2V ^ dx < ST\Xf x \ 2 . 

e^O 0<t<T Jo JdD 



From Proposition 12.41 item 3, we can choose A small enough so as to make the latter term 

arbitrarily small. So we complete the proof. □ 

Proof of Theorem HOI 1) From [fTT|). we only have to check that ([31]) holds under P e *. This 
follows from Theorem 12.71 and the estimate (obtained with Lemma 12 - 1[) 



lim E f 



SU P I / (fe - f)( T Xt/eU)dr\ 
}<t<T Jo 



'<T\f £ -f\ 



-0<t<T Jo 

The same argument holds for (|33|) . □ 
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E Proofs of subsection 12.51 

Proof of Provosition \2.1(A The statement (f3l)f is quite classical. The reader is referred to [T^l Ch. 
2] for an insight of the method and to [T9J, Prop. 4.3] for a proof in a more general context. □ 
Proof of Provosition \2.1l\ In what follows, for each i = 1, • • • , d, (<p l „) n stands for a sequence in 
C such that D<p l n -> C in L 2 (Q) d asn-> +00. 

Let us first focus on (f39|) . Fix X £ R d whose entries are denoted by (Xi)i<i<d- We have: 
D{X i ip i n ) = XiDtpi -> Xid = CX in L 2 (Q) d as n -> +00 and 

X*AX=M[(X + CX)*a(^ + C*)l = lim M[(Jf + DpQ<^))*a(^ + 

n— >-t-oo 

> inf M[(X + Dip)*a(X + Dip)}. 
Conversely, from Lemma IE. II below, we have: 

(61) VY £R d , M[{Y + CY)*aCX}= lim M[(Y + CY)*aD(Xi<pi)] = 0. 

n — > + oo 

The above relation and Lemma [Ej] again yield M[(X + CX)*a(Dip — CX)} = for any tp £ C. 
So, for every ip £ C, we have: 

U[(X + Dip)*a(X + Dip)} =M[(X + CX + Dip - CX)*a(X + CX + Dip - CX)} 

=M[(X + CX)*a{X + CX)} + 2M[(X + CX)*a{Dip - CX)} 

+ M[(Dip - CX)*a(Dip - CX)} 
>M[(X + CX)*a(X + CX)} 

so that (J2U) follows. By the way, (fBTj) proves that A also matches M[(I + C*) a }- 

Now we prove AI < A. Fix X £ R d . From © and Cauchy-Schwarz's inequality, we get 

X*AX = M[(X + CX)*a(X + CX)} > AM[\X + CX\ 2 } > A\M[X + CX} | 2 = A\X\ 2 , 

since M[CX} = 0. The estimate AI < A follows. 

Now we prove that T = M[(I + C*)l] coincides with the orthogonal projection Mi [(I + C*)7]- 
Remind that *y can be rewritten as 7 = ae\. So we just have to establish the relation 

(62) Mi[(I + C>ei] = M[(I + C*)oei]. 

Proof of (|6"2"| . Because of the ergodicity of the measure fj, (2. of Definition 1 we stress that a 
function tp £ L 2 (fl, Q* , fj,) invariant under the translations {t x ; x £ R x {0} d_1 } must be constant 
and therefore satisfies Mi[i/i] = M[-0]. So we just have to prove that the entries Wli[(ei + C e i)* ae i] 
are invariant under the translations {t x ; x £M.x {0} d_1 }. To that purpose, we only need to check 
that 

M[M 1 [(e 4 + Ce 4 )*aei] J Di V ] =0 
for any i = 1, . . . , d and ip £ C. By using Lemma |E. 21 ii below, we get: 

M[M![( ei + Ce i Tae 1 ]D 1 ip} = M[(e, + Cei)* ae l M 1 [D l ip\) = M[(e< + Ce i )*ae 1 D 1 M 1 [ V >}} . 

Since L> fc Mi[<p] = for k = 2, . . . , d (see Lemma EH i), we have eiDjMi[tp] = DMi [<p]. Wc 
deduce 

M[M a [(e 4 + Ce i )*ae l }D 1 ip\ = M[(e, + Ce i )*aDM 1 [if}} . 

Since Mi[y>] £ C (Lemma IE. 2l ii). the latter quantity is equal to (Lemma lE.l[) and we complete 
the proof. Note that the above computations also prove: Ti — M[(ei + C e i)* ae i] = An > A. □ 

Lemma E.l. The following relation holds: 

(63) VX £ R d , V^GH, M[(X + CX)*aDi(>} = 0. 
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Proof. Since bi — hD k ai k , the weak form of the resolvent equation (fH)|) associated to / = bi 
reads, for any tp E M: 

X(u\,ip) 2 + (l/2)(a jk D jU {,D k xp) 2 = (l/2)(D fe a ife , ^) 2 = -(l/2)(a ik ,D k ip) 2 . 

By letting A go to and by using we obtain: (l/2)(a jk Cj, D k~4>) 2 = -{l/2)(a ik ,D k ip) 2 . 

We deduce M[(S tj + C)) 

aj k D k xjj\ = 0, which means nothing but 
(64) M[( ei + Cei)*al?^]=0. 

The result follows by linearity. □ 

Lemma E.2. The projection operator Mi saisfies the following elementary properties: 
i)Vk = 2,...,d and G Dom(D fe ), D fe Mi[y>] = Mi[D k ip] = 0, 
n) Vcp e C, Mi[(p] e C and Mi[Di<p] = DiM^ip], 
in) yk = 2,...,d andyip.ip e Dom(D fc ), M^Dkipip] = -M^ipDkip]. 

Proof. The properties i) and ii) are easily derived from the identities MifT^ip] = Mi[y>] for any 
x e {0} x R d ~\ T X M 1 = MiT x for anyielx {0} d ~\ and Mi [tp * p] = M^tp] * p for any 
ip € L°°(ft) and p € C^°(M d ). iii) results from i). Details are left to the reader. □ 

F J-topology 

We summarized below the main properties of the Jakubowski topology (J-topology) on the space 
D([0, T];R) (set of functions that are right-continuous with left-limits on [0, T]) and refer the 
reader to [TT] for further details and proofs. We denote by V the set of functions v : [0, T] — > R 
with bounded variations. The J-topology is a sequential topology defined by 

Definition F.l. A sequence (x n ) n in D([0,T];R) converges to xq S D([0, T];R) if for every 
e > 0, one can find elements (v n .e)neN C V such that 

1) for every n € N, supr T i \x n — v n ,e\ < £, 

V/ : [0,T] — > K continuous, f Q f(r)dv„ ie (r) — > J Q T f(r)dvo. £ (r) as n — > +oo. 

By gathering [UJ Th. 3.8] and [HI Th. 3.10], one can state: 

Theorem F.2. Let (V a ) a C £)([0,T];R) be a family of nondecreasing stochastic processes. Sup- 
pose that the family (V a (T)) a is tight. Then the family (V a ) a is tight for the J-topology. More- 
over, there exists a sequence (V n ) n C (V a ) a , a, nondecreasing right- continuous process Vo and a 
countable subset C C [0, T[ such that for all finite sequence {t\, . . . ,t p ) C [0, T] \ C, the family 
{V n (tx), . . . , V n (tp)) n converges in law towards (Vo(ti), ■ ■ ■ , Va(t p )) n . 

Equip the set V+([0, T]; R) of continuous nondecreasing functions on [0, T] with the J-topology 
and C([0, T];R) with the sup-norm topology. We claim: 

Lemma F.3. Let (V n ) n be a sequence in converging for the J-topology towards Vq € V+. 
Then (V n ) n converges towards Vq for the sup-norm topology. 

Proof. This results from Corollary 2.9 in [TT] and the Dini theorem. □ 
Lemma F.4. The following mapping is continuous 

(i,v)eC([0,r];I)x¥+([0,r];l)H f x r dv{r) e C([0,T];M). 

Jo 

Proof. This results from Lemma IF.3I and the continuity of the mapping 

(i,«)eC([0 ) r];R)x¥+([0,I1;l)H f x r dv{r) e C([0, T}; R), 

Jo 

where both C([0,T];R) and V+([0,T];R) are equipped with the sup-norm topology. The reader 
may find a proof of the continuity of the above mapping in the proof of Lemma 3.3 in [17] (remark 
that, in [TT], the S-topology coincides on C([0,T];R) with the sup-norm topology). □ 
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G Proof of the tightness (Proposition 12.121) 

We now investigate the tightness of the process X £ (and K £ ). Roughly speaking, our proof is 
inspired by [16( Chap. 3] and is based on the Garsia-Rodemich-Rumsey inequality: 

Proposition G.l. (Garsia-Rodemich-Rumsey's inequality). Let p and "J be strictly in- 
creasing continuous functions on [0, +oo[ satisfying p(0) — 'J'(O) = and lim^oo ^S(t) = +oo. 
For given T > and f E C([0, T]:R. d ), suppose that there exists a finite B such that; 



(65) r f T H i9it t 9 ^) 

Jo Jo V p(\t-s\) J 



dsdt < B < oo . 



Then, for all0<s<t< T: \g(t) - g(s)\ < 8 f* s ^{AB/v?) dp{u). 

To apply Proposition [GUI it is necessary to establish exponential bounds for the drift of X e . 
Indeed, suppose that we can prove the following exponential bound: for every < s, t < T 
(66) 

W* exp(n\J [^b j -d Xi V(X £ )a ij ](T X e/ e uj)]dr+ J a l:j {T X ,/ £ uj)dK £ \^ < 2 exp (Cn 2 (t - s)) . 

for some constant C > depending only on A (defined in ©). Then we can apply Proposition 
[Gil as detailed in 16, Ch. 3, Th 3.5] (set p(t) = s/i, ip(t) = e* - 1 and ^A" 1 ^) = ln(i + 1) in 
Proposition [GT} to obtain 

Proposition G.2. We have the following estimate of the modulus of continuity 
(67) 

E e *( sup \f [-b j -d Xi V(X*)a ij ](Txe/sU)]dr+ [ aij (T X * /e u) dK e r \) < C^ln^ 1 ), 

v \t-s\<S;0<s,t<T Js e Js ' 

for some constant C that only depends on T, A. 

We easily deduce the proof of Proposition 12.121 we first work under P e * . Let us investigate 
the tightness of X e . Observe that 

x t' e = x j+ [- b 3( T x-/eu)-d Xi V(X £ )a lj (T X . /e uj)]dr+ aij(Txf/ e w) dK*+ / a ji (T X e./ e w)dB*' 1 
Jo £ Jo Jo 

The tightness of the martingale part follows from the boundedness of er and the Kolmogorov 
criterion. The tightness of the remaining terms results from Proposition IG.2I So X e is tight 
under P e *. 

Let us now investigate the tightness of the family (K e ) e . From Lemma |2"7T1 we have W* [Kj] = 
T J gD e ~ 2V{ - x) dx. Theorem E2 ensures that {K e ) £ is tight in D([0,T];R+) (remind that K £ is 
increasing) . 

To sum up, under P", the family (X £ ,K £ ) £ is tight in C([0,T];D) x D([0,T];R+) equipped 
with the product topology. From ([IT]), the family is tight in C([0,T];D) x D([0,T];R + ) under 
P e . □ 

We have thus shown that the proof of Proposition (|2.12|1 boils down to establishing (|66|) . So 
we now focus on the proof of (f66|) . We want to adapt the arguments of [16, Chap. 3]. However, 
the situation is more complicated due to the pushing of the local time when X e is located on the 
boundary dD. Our idea is to eliminate the boundary effects by considering first a truncated drift 
vanishing near the boundary: fix lu E £1 and a smooth function p E C£°(D) satisfying p(x) = 
whenever x\ < 9 for some 9 > 0. For any e > and j = 1, . . . , d, define the "truncated" drift 

(68) KA X ^) = e —^d Xi {e- 2V ^a ij {r x/ ^)p{x)), 
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which belongs to C£°(D). Our strategy is the following: we derive exponential bounds for the 
process J Q b p j(X^,w) dr. These estimates will depend on p. Then we shall prove that we can 
choose an appropriate sequence (p n ) n C C£°(D) preserving the exponential bounds and such that 
the sequence ( f Q b £ Pn j(X^, uj) dr) converges as n — ► oo towards the process involved in (f66|) . 

The exponential bounds are derived from a proper spectral gap of the operator C v ■ +nb s p •• 
with boundary condition 'Yi(T x / E u))d Xi ■ — on dD. The particular truncation we choose in (|58|) 
is fundamental to establish such a spectral gap because it preserves the " divergence structure" of 
the problem. Any other (and maybe more natural) truncation fails to have satisfactory spectral 
properties. 

So we define the set 

C 2 * = {/ G CZ(D) ni (T x/e u)d x J(x) = for x G dD} 

and consider the Hilbert space L 2 (D; e~ 2V ^dx) equipped with its norm \ -\rj and its inner product 
(-, -) D . Given k > and oj G Q, let V£ K G C°°([0, T] x D) n C b 1,2 be the unique solution of 

d t rj K = C £ v^ K + K KA^ K + !) 011 [°. T l x D ' 7^/^)3*^'" - on [0, T] x dD 



with initial condition ^'"(O, •) = on D (see Lemma [A.l[l . Then uJ,' K = ip^J K + 1 € C h ' is a 
bounded classical solution of the problem 

(69) d t u 6 j K = C v u e j K + Kb pJ u e j K on [0, T] x D, 1^,^)8^ = on [0, T] x 9£>, 

with initial condition m £,k (0, •) = 1 on D. Lemma \A . 2 1 and a straightforward calculation provide 
the probabilistic representation 



u^{t,x)=W x * / K &^.(A7, w )exp( / K6^(X«,a;)du)dr 



exp(«; / b pj {X e rl uj)dr) 



Lemma G.3. For each uj G f2, we /icrae i/ie estimate \ufj K (t, -)\ 2 D < e 2t7T ^'" (0 < t < T), where 
ir^j K = sup(^, £y</> + Kb e p j(j))D o,nd the sup is taken over {(f) G C^ ,e , = 1}. 

Proof. We have: 

9t|ur(«,-)lD=2«' K ,9t<' K (t,-))D 

=2«;«,^<;' t + ^^(t, 0)d < 27C«K»(t, Oil,. 

Since |m^' k (0, ■)!!, = 1, we complete the proof with the Gronwall lemma. □ 
Proposition G.4. For any k > 0, e > and < s, t < T 



' [ exp ( | k / b e p<j (X*,uj)dr\)] < 2exp(CK 2 (i-s)), 



/or some constant C that only depends on A and sup xe jj |p(x)|. 

Proof. By stationarity (resulting from Lemma l2.1[) and Lemma lG.3[ we have 

E s *[exp(«; J b p] {X £ rl uj)dr)} <E £ *[exp( K J b pJ (X e r ,uj)dr)] 
(70) =M* D [u e J K (t-s,x)] <M\u s J K (t-s,-)\ D < M[exp((i - s)<'")]. 

It remains to estimate 7rf/ K . For any function <fi G Cy £ such that \<j>\jj = 1, we have 

(6p J -(-,cj),0 2 ) D = -(aij(T./ e w)p(l>,d 3 . i (f>) D < A -1 sup ^(x)!!^^^ = Cl^^Ji 
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where we have set C = A 1 snp xe Q \p(x)\. As a consequence (the sup below are taken over 
{0 6^,1012, = !}) 



TT £ J K = SUp((f>, £ £ y(f> + Kb S pJ (j>) D 

< sup{ - {l/2)(a lj (T. /e uj)d Xz <t> 1 d X] (j)) D + n(b e pj {-, w), cj) 2 ) D } 
(71) <sap{-(A/2)\B x 4>\% + KC\B x 4>\ D } <k 2 C 2 /(2A). 

The last inequality is obtained by optimizing the expression — (A/2)x 2 + kCx with respect to the 
parameter x £ R. Gathering (|70|) and (j7"Tj) then yields 

E e *[exp(K^ 6^(X*,w)<2r)] < exp (C"K 2 (t - s)) 

where C — sup^g^ |p(a;)| 2 /(2A 3 ). We complete the proof by repeating the argument for —b e p j 
and using the inequality exp(|x|) < exp(— x) + exp(x). □ 

As explained above, we can replace p in Proposition lG.4l with an appropriate sequence (p n )n C 
C£°(D) so as to make the sequence ( J Q b pn j(X £ ,uj)dr) n converging as n — > oo towards the 
process involved in (JHH). Let us construct such a sequence. For each n £ N* , let us consider the 
piecewise afhne function p n : D — > M defined by: 

1 1.12 

Pn(%) = if Xi < — , p n {x) = n(x\ ) if — < x\ < — , and 1 otherwise. 

n n n n 

Note that p n is continuous and sup^gf) |p„(a;)| < 1. With the help of a regularization procedure 
and Lemma I2.1[ one can prove that Proposition IG.4I remains valid for p n instead of p, where 



(72) 



* 1, 



b pn!j (X e r ,u)dr = / [-bj(Tx. /e u) - d Xi V(X £ ) aij (T X?/s io)]p n (X £ )dr 
Jo £ 

+ / o,ij(T X s / e w)nlM.2,(I r e ) dr. 
Jo 



You can obtain the latter expression by expanding (|68j) with respect to the operator d Xi ■ 
Since sup^g^, |/5n(-^)| = 1 for each n, we deduce 

(73) Vn e N,V0 < s,i < T, E £ *[exp(|«;^ b pn j(X e , w) dr|)] < 2 exp (Cn 2 {t — s)) 

for some constant C only depending on A. Now it remains to pass to the limit as n — > oo in (|73p . 
From [U Prop 1.19], J Q aij(Tx^/ E cj)nTl^i.2^(X £ )dr converges a.s. towards J Q aij(T X e/ e u)dK e as 
n — * oo. Fatou's Lemma (as n — > +oo) in ([73^) then yields ([5o]) . So we complete the proof. □ 
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